SECuRtTy  cu  ^SSt  riC  A'T’iOn  TmiS  »»AOe  '**Ti#n  0«f*  Entmfa. 


REPORT  DOCUMENTATION  PAGE  j 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  report  number  2  GOVT  ACCESSION  N^ 

NRL  Memorandum  Report  4564  i 

)  RECIPIENT'S  CaTacOG  number 

4  TiT<.E 

NONLINEAR  ST.ABILIZATION  OF  THE  FARLEY-  j 

BUNEMAN  INSTABILITY  BY  STRONG  E  X  B  | 

TURBULENCE* 

5  Type  of  report  a  period  covcrco 

Interim  report  on  a  continuing 
NRL  problem. 

a.  performing  drc.  report  number 

7  AUTHOR/' *>» 

8.  contract  or  grant  NuMBCRr*; 

M.  J.  Keskinen 

9.  RERFoRMInO  CRGAni  2  A'^iON  NAmE  AnQ  ADDRESS  j 

Naval  Research  Laboratory 

Washington,  DC  20375 

1  10  program  element,  project,  task 
AREA  a  work  unit  numbers 

61153N;  RR033-02-44; 
47-0883-0-1 

n  control;-ing  oFPice  name  and  address  I 

Office  of  Naval  Research  i 

12.  REPORT  DATE 

July  31,  1981 

Arlington,  VA  22217 

1  1].  number  of  pages 

i . — 

UNCLASSIFIED 


is«  declassification  downgrading 
schedule 


Accession  For 


I*.  OlST’WiSu  ^»ON  ST^ATEmEnT  rot  thtM  Report f 

Approved  for  public  release;  distribution  unlimited. 


NTIS  CRAftI 
DTIC  TAB 

Unannounced  □ 
Justification _ 


!7  OISTNiBu'^'On  Statement  at  '/I#  *ntmf4  \n  Sfocli  20.  it  dliimtmfxt  /row 


By - 

Distribution/ 


Dist 


Availability  Codes 
]Avnil  niid/or  . 
Special 


'i.  SuRWL  En£n  tary  nO*ES 


fi 


♦This  work  was  completed  while  the  author  was  an  invited  guest  at  the  Centre  de  Physique 
Theorique,  Ecole  Polytechnique,  Palaiseau,  France. 


1$  KEY  WORDS  '  Confinu*  on  ltd*  it  n«c««t«rv  mnd  tdmntlfy  Sy  block  numbmr) 

Strong  turbulence 
Equatorial  electrojet 
Nonlinear  saturation 
Farely-Buneman  instability 


20  abstract  'Cofittnum  on  rovrao  tfd*  //  n«ee«««0'  wnrf  Idmntlty  by  block  nuatbor) 

~  It  is  shown  that  through  nonlinear  mode  coupling  processes  long  wavelength  low 
frequency  strong  E  X  B  turbulence  can  stabilize  short  wavelength  high  frequency  Farley- 
Buneman  modes  in  a  weakly  ionized  low  pressure  convecting  plasma.  Favorable 
comparisons  are  made  with  experimental  observations.  *- 


DO , 1473  EDITION  OF  1  NOV  as  IS  OBSOLETE 

S  N  uo  2-  3  l«-  »60  I 


securitt  Classification  of  this  race  iWh»#i  o«<« 


1 


NONUNEAR  STABILIZATION  OF  THE  FARLEY-BUNEMAN 
INSTABILITY  BY  STRONG  E  X  B  TURBULENCE 

It  is  well  known  that  in  the  absence  of  a  magnetic  field  the  two  stream 
instability  can  occur  in  a  homogeneous  plasma  when  the  electron  drift  velo¬ 
city  with  respect  to  the  ions  exceeds  the  electron  thermal  velocity.^ 

Farley^  and  Buneman^  have  shown  that,  in  the  presence  of  a  magnetic  field,' 
the  electron  drift  velocity  with  respect  to  the  ions  has  only  to  exceed  the 
ion  acoustic  velocity  to  generate  unstable  waves  traveling  perpendicular 
to  the  magnetic  field.  We  will  consider  the  nonlinear  evolution  of  the 
Farley-Buneman  instability  in  a  low  0,  weakly  ionized,  convecting  plasma 
which  is  subjected  to  a  magnetic  field  Bx,  an  electric  field  E^z  and  a 
density  gradient  (Bn^/dz)z.  Differences  in  the  collision  frequencies 
(Vi/Q^  ^  1,  I)  of  the  ions  and  electrons  with  the  background  neu¬ 

tral  gas  results  in  the  formation  of  a  cross  field  current  'd^y  from  the 

^  £  X  electron  drift.  For  weak  currents  long  wavelength  field 

o 

aligned  fluctuations  in  density  6n  expCi{ky-ajt)]  have  been  found  by 

Simon^  and  Hoh^  to  be  linearly  unstable  when  E  'Vn  >0.  In  the  nonlinear 

o 

regime  this  ^  x  B  gradient  drift  instability  evolves  into  an  isotropic  two- 
dimensional  strongly  turbulent  state  in  the  plane  perpendicular  to  the  mag¬ 
netic  field  as  shown  previously.^  For  stronger  currents  such  that  V  >  C 

Q  S 

the  Farley-Buneman  instability  will  develop  at  shorter  wavelengths.  These 
long  and  short  wavelength  modes  can  coexist  simultaneously  (see  Fig.  1)  with 
the  former  usually  occurring  before  the  latter.  Previous  studies  of  the 
nonlinear  evolution  and  saturation  of  these  short  wavelength  Farley- 
Buneman  modes  have  invoked  quasilinear  effects,^  resonance  broadening,^  and 

mode  coupling.^  These  works  have  neglected  the  effects  of  the  strong  large 
scale  background  E  x  ^  turbulence.  In  this  Letter  we  show  that  the  long 

wavelength  ^  x  ^  turbulence  can  stabilize  the  short  wavelength  high  fre¬ 
quency  Farley-Buneman  instability.  Although  the  following  discussion  is 
Manuscript  submitted  May  12,  1981. 
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applicable  to  any  weakly  ionized  low  3  current  carrying  plasma  convecting  in 
regions  of  ^  x  ^  turbulence  it  has  direct  bearing  on  density  irregularities 
in  the  equatorial  electro jet  ionospheric  plasma. 

The  basic  equations  for  the  electron  (N^)  and  ion  (N^,)  fluids  in  a  low 
3,  werkly  ionized,  collisional  plasma  can  be  written  <' 

oN/dt  +  7»NV  »  0  (1) 

e.'JCE  +  V  X  B)  +  IVN  +  Nm  v  V  =  0  (2) 

3N/3t  +  V.w  =  0  (3) 

~i 

(5/dt  +  V  *7)^111.7,  =  eN(E  +  V.  x  B)  -  T7N  -  Nm.v.V.  (4) 

(5) 

where  we  have  assumed  quasineutrality  (N^  =“  M) ,  isothermality  (T^  =*•  T^) , 

electrostatic  fluctuations  and  neglected  electron  inertia. 

Linearizing  equations  (l)-(5)  with  N  =  n^  +  n,  etc.  and  assuming  fluctua¬ 
tions  of  the  form  n,  9,  v^,  «  exp|i£k,^y  +  k^z  -  +•  iY(k)t]}with 

k*^  ■  0  and  k^L  »  1  we  find®  for  the  frequency  and  growth  rate  in  the  ion 
frame 

•^kr  “ 

''(k)  -  C’i'/ (1  +  ’I*)]  I  (V^/L)cos20  +  ^ 

where  ^  ^  “  (1/n^)  (Sn^/iz) ,  «  2kT/m^  and  0  is  the  angle 

defined  by  Jk  and  From  the  expression  for  the  growth  rate  Y(k)  in 


2 


(7),  we  note  that  at  low  frequencies  (long  wavelengths  and  weak  currents) 

such  that  tL'  /v  <  (Q^/v^)  (1/kL)  the  E  x  B  gradient  drift  term  will  dominate 

with  all  modes  with  k  <  k^  «  Q  (0^/v^) (V^/L) (v^/c^^) cos0  unstable.  At 

higher  frequencies  (short  wavelengths  and  strong  currents)  Farley-Buneman 

modes  will  become  unstable  if  but  with  no  critical  wavelength.  In  / 

the  equatorial  electrojet  plasma  V,  (t)  is  time  dependent  varying  from 

d 

V  <  c  to  V,  ^  c  over  a  time  interval  At  »  Y(k)  As  a  result  the 
Farley-Buneraan  instability  will  be  excited  in  strong  ^  x  ^  turbulence.  Since 
these  waves  are  nondispersive  they  will  interact  strongly.  However  as  their 
amplitude  increases,  they  do  not  steepen  appreciably^  but  are  unstable  to 
perturbations  perpendicular  to  their  propagation.  It  must  be  noted  that  this 
fluid  approximation  is  valid  for  ui  ^  v^.  At  higher  frequencies  kinetic 
effects,  e.g.,  ion  Landau  damping,  will  become  important  and  introduce  a  high 
frequency  cutoff.  For  lower  frequencies  both  fluid  and  kinetic  treatments 
are  identical.^ 

By  writing  N(x,t)  =  n  +  L  n(k,aj)expC-i(k*x-u)t)  ] ,  $(x,t)  =  cp  +  Z  cp(k,a)) 
^  o  kuj  ^  ^  ^  ku) 

expC-i(k  •x-u)t)],  etc.,  expanding  eqs.  (l)-(5)  in  the  small  parameter 

^  ^  ^  considering  high  frequencies  such  that 

uu^/v^  >  (C^/v^)  (1/kL)  ^  10“^  (for  \  ^  3m,  L  ^  6  km,  *  10^)  we  find  to 

second  order  in  n^  (n/n  «  1) 

D(k,u))ri (k,ou)  »  /d^k' duu' V(k,k'  )  n(k'  ,u)'  )  n(k-k'  .oj-nju'  )  (8) 

where 

D(k,ii))  *  a)-k*V,(l  +  ’!;)  •-ivw^(l  +  ij^)  v  +  ii);k^c  ^/v. 

i  Si 

is  the  Farley-Buneman  dielectric  and 
V(k,k'  ,«ju,uj'  )  »  -(x  X  k'  ‘k/k'  2)  [(i  +  (v^/a^)k'  -V^  +  k'  ‘x  x  V^+ik' 
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By  neglecting  the  nonlinear  term  on  the  right  hand  side  of  eq.  (8)  we 
recover  the  linear  result  from  D(jc,uu^^  +  iyCk))  ■  0  giving 

A/ 

and  y(k)  »  (^/l  +  ’«^)(l/v^)  -  k^c^^] . 

solve  eq.  (8)  for  the  high  frequency  short  wavelength  component  of 
n(k,a>)  by  considering  its  mode  coupling  co  the  low  frequency  long  wavelengtft 
well  developed  strong  ^  x  ^  turbulence.  Let  a  Farley-Buneman  wave  be  denoted 
by  (kj,a)j)  and  a  turbulent  ^  x  ^  mode  by  Physically,  when  a 

Farley-Buneman  mode  (k  ,uj-)  grows  to  such  a  level  that  it  can  couple  with  the 

X. 

E  X  B  turbulence  (k__,'jj-_),  a  beat  wave  component  (k^.  i  will 

^  ^11  II  i  II 

appear  which  in  turn  can  beat  with  (k-_,a)__)  to  affect  (k.jU)^).  The  evolu- 

I  II  *^1  I 

tion  of  nCk-.uJ.)  can  then  be  written 
I 


D^(k-,uu-)n^(k^,iu-)  =  /d^k' duu' V(k, ,k'  ,0)-,^  )n^^(k'  )n^  ^^(k^-k'  ) 

^I  I  />^X  I  j  ^I  ^  '  I  ^I  I 


(9) 


The  beat  wave  (k.-k-- evolves  according  to: 
^1  I  II 


D^"^^(k^-k,uj^-uj)n^“^^(k_-k,u)^-u))  =  / d^k'  dw  VCk^-k.k'.uu^-iJu ,uo'  )n^^(-k'  ,-uu'  ) 

0>»l  ^  I  ^I  I  J  ^  ^  L 


n  (k 

^1  I 


(10) 


where  we  have  used  V(k,k^  “  V(k,-k^  Substituting  eq.  (10)  into 

eq.  (9)  we  find  the  nonlinear  dispersion  relation  to  lowest  order 


D^(k^,uJj)n^(kj,a;j)  -  0 


(11) 


T  T 

where  D  (k-,d)«)  *  D  (k_,a)-)  +  6D  (k-,jj-).  The  nonlinear  part  of  eq.  (11) 
I  I  I 

can  be  written 


/•  V(k  ,k'  ,uj  ,u))V(k^-k'  ,k 

6D^(k_,a)T)  -  -  /d^k'du)'  - - - - - - - -  I^\k'  ,ou' ) 

^  D  (k--k 

^  I 


(lla) 
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I 


b 

i 

( 


» 


ki 


■  i 


where  *  <  n^^(k,a))n^^(-k,-u)) >  is  the  power  spectrum  of  the 

£  X  ^  turbulence.  It  has  been  previously  shown^^  that  I^^(k,ui))  can  be  calcu¬ 
lated  using  the  direct-interaction  approximation  of  Kraichnan^^ 


lu)-u)(k)  +  r^^(k,aj)  1  -  j  J'd^^  dw  lw(k,k'  )l^I^^(k'  ,u)'  )I^^(k-k'  ,ujW  ) 

(12) 

TT  r  w(k,k-k'  )w(k-k'  ,k)I^^(k'  y  ) 

r'-^Ck.uu)  «  -  Jd\'dw  - — -  (13) 

uu-uu'  -U)(k-k'  )+r^^(k-k'  ,uuV  ) 


with  w(k,k^  )  »  VCkjk^  )  +  V(k,k-k^  )  and  aj(k)  »  o),  +  iY(k)  where 

yOc)  =  (i///l  +  (V^/L)cos^0  -  jc^Cg^/v^],  In  eq.  (12)  r^^(k,cu)  is  the 

self-damping  of  the  long  wavelength  E  x  B  fluctuations  (k-_,u)__).  We  now  pro- 

/N*  /W  11 

ceed  to  solve  equations  (11)-(13). 

First,  we  note  that  since  the  interacting  waves  considered  here  are  non- 

dispersive  (see  eq.  (6))  the  beat  wave  dielectric  D^"^^(k--k-^ 0  and 

I  II 

the  right  hand  side  of  eq.  (11a)  diverges.  As  a  result  we  replace 

^^(k_-k^  )  by  its  renormalized  value  ^^(kT.-k^  )  + 

6D^  ^^(k^-k^  6D^  order  to  study  the 

nonlinear  saturation  we  solve  eq.  (11a)  for  Im  6D^  =  -ir(k_  ,uj- (k.) ) .  From 

^1  1  ^1 

previous  studies'* ’  the  steady  state  solution  of  eqs.  (12)-(13)  can  be 
vrritten 


I^^(k,uu)  »  l^^(k)  (2TT)”^(r^^(k))“^exp[-(uu-u)(k))2/2(r^^(k))2]  (14) 

where  u)(k)  =  k*^/(l+'i/),  r^^(k)  >■  3.4  n  ^(v  /Q  )k^V,  (I^^(k))  ,  and 

I^^(k)  *  Ik"^  is  isotropic^ » ^ ^  with  n  3-4  while^^  I  «=  V,^,  m  ^  2.  Substitu- 

d 

ting  eq.  (14)  into  (11a)  and  assuming  that  k^  ^  k^^  <  k^  s  k,  uu'  < 

j)-  *  0),  ^^(k--k'  ,uj  -0)'  )  ^  D^(k  ,0)-)  +  0(k'  /k-)  ^  6D^(k_,a)-)  we  find 

1  ^1  I  <^i  1  1  ^I  I 


EQ 
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[r(k)]2  »  yd2k'v(k,k'  )V(k-k'  y  )I^^(k'  ) 

^  X  k'  •k)2  (k'  •^)^lk'  j-**  I^^(k'  ) 

“  (v^/n^)2k2v^2  yde' k' dk'  Sin2(e-e'  )cos2e'  I^^(k'  )  (15) 

where  r(k)  >>y(^)  has  been  assumed  and  0^  and  9  are  Che  angles  made  by  k''  and 
k,  respectively,  with  V,.  In  evaluating  the  quantity  V(k,k^)  V  (k-k^,k) 
in  eq.  (15)  we  have  kept  only  terms  proportional  to  where 

^i^^i  Experimental  studies^^  as  shown  in  Figure  2(a)  indicate  that  r(k) 

is  approximately  independent  of  angle  9.  This  allows  the  replacement  of  r(k) 
in  eq.  (15)  by  its  angle  averaged  result  (2tt)  ^^d0r(k)=r(k)  giving 

r(k)  =  (v^/n.)(y)  kV^  <  ln/n^l2^^  (16) 

where  <  \nln^\^  dk  I^^(k)  =  2Tr  Jdk  k  I^^(k).  Using  22, 

2Tr/k  ^  5m,  =  4  x  102m/sec,  ^  3.6  x  102m/sec,  and  <  [n/n^P^^  >'’^  =“  0.01, 

eq.  (16)  gives  for  the  spectral  width  r(k)/2T=**  10  Hz  (cf.  Fig.  1) .  Figure  2(b) 
shows  that  the  scaling  of  r(k)  «  k  from  eq .  (16)  is  in  reasonable  agreement 
with  experimental  results"- which  indicate  that  the  spectral  broadening 
iiu,^  «  r  (k)  oc  k^  ■  . 

The  nonlinear  dispersion  relation  in  eq.  (11)  for  the  short 
wavelength  Farley-Buneman  instability  in  the  long  wavelength  ^  x  ^  turbulent 
background  can  then  be  written  explicitly 

1  h 

+  vp)  '■v^^+  i!i/k2c^2/^^^^  .  0. 

(17) 

Separating  uu  ■  +  iY(k)  we  find 

1  ^ 

Y(k)  *  +  il»)]  I  (k*j^)^/v^  -  k2cg2/v^[  -  •j(v^/r.^)kV^<l  n/n^l  2j^> 

6 


(17a) 


For  complete  stabilization  of  the  fastest  growing  linear  mode  (8^0)  the  long 

wavelength  background  fluctuation  level 

<  In/n  [2  ^  2  V  Q  -‘(1  +  /v.  -  k^c  2/v  kV  ). 

o^^  ee  ni  sru 

with  kV,/v.  ^  0.1,  k^c  ^/v.kV,  "^0.05  we  find  n/n  ^  0.002  which  is  consis- 
a  1  s  1  a  o 

tent  with  available  experimental  estimates  . 

i 

Part  of  this  work  was  completed  at  the  Centre  de  Physique  Theorique, 
Ecole  Poly technique,  Palaiseau,  France.  We  wish  to  thank  the  Centre  National 
de  la  Recherche  Scientifique  for  their  invitation  to  visit  the  Ecole 
Poly  technique.  Useful  discussions  with  R.  Pel!)  at  of  the  Ecole  Poly  technique 
and  M.  Crochet  and  C.  Hanuise  of  the  Universite^  de  Toulon  et  du  Var  are 
acknowledged. 
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Doppler  Shift  (Hz) 

Fig.  1  —  Simultaneous  equatorial  electrojet  ionospheric  plasma  density  fluctuation  power 
spectra  (vertical  axis  I(k,cj)  vs.  frequency  w  in  Hz  (horizontal  axis)  from  Ref.  14  for  several 
radar  backscatter  observation  frequencies;  the  radar  frequencies  29.0  MHz,  9.0  MHz, 
6.6  MHz,  5.6  MHz  correspond  to  wavelengths  5.1m,  10.5m,  16.6m,  22.7m,  and  26.7m, 
respectively. 


(a) 


V  5.2  m 
□  10.8  m 

0  21.7  m 
A  39.3  m 


^  i 

A  A 
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6 

A/? 


0.5 


Fig.  2  —  (a)  Power  spectrum  width  vs.  wave  vector  angle  made  by  k  and  the  electrojet 
electron  drift  velocity,  at  several  wavelengths  from  Ref.  14;  note  approximate  isotropy  in 
angle,  (b)  Power  spectrum  width  vs.  wave  number  k  at  d  -  15°  from  Ref.  14;  solid  line  is 
best  fit  a  k  . 
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